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5.0 MODEL VERIFICATION 


In order to establish the validity of the EMA model, comparisons 
are made between results obtained from lab tests on the actual EMA and 
those from equivalent computer simulations. The computer program has 
several run options at. described earlier in Section 4.5, for this 
reason, the following verification schedule will be followed: 

1. Verification of the Power Conditioner and Machine Model 

• Motoring 

• Regeneration 

• Plugging 

2. Verification of the Fourteenth Order EMA Model 

3. Verification of the Tenth Order EMA Model 

Table (5.0-1) catalogs the various computer runs which are used in 
this verification procedure. The results of these runs are presented 
in the twelve appendices: A.1-A.3, B.l, C.1-C.3, D.l, E.1-E.4. In 

terms of the aforementioned verification program, these appendices 
are organized as follows: 

1. Power Conditioner and Machine Model 

• Motoring: Appendices A.l, A. 2, and A. 3 

• Regeneration: Appendix B.l 

• Plugging: Appendices C.l, C.2, and C.3 

2. Fourteenth Order EMA Model: Appendix D.l 

3. Tenth Order EMA Model: Appendices E.l, E.2, E.3, and E.4. 

The tests on the actual hardware were performed by Delco Electronics 
and are described in References (2) , (3) , and (4) . 

A summary of each computer run is given on the first page of 
each of the twelve appendices. The results of these computer simu- 


TABLE 5.0-1 
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lations are plotted by a CalComp pen type plotter. Tables 5.0-2 and 
5.0-3 contain listings and descriptions of the CalComp plots. Figures 
(A. 1-1), (A. 2-1), (A. 3-1), (B.l-1) , (C.l-1), (C .2-1) , (C.3-1), (D.l-1) 
(E.l-1) , (E .2-1) , (E.3-1) , and (E.4-1) give the status of the logic 
signals versus time for each of the twelve computer runs. These (5.0-4) 
signals are defined in Table (5.0-4). With these logical 
variables the reader can easilty follow the status of each solid state 
power conditioner element as time progresses. 

A few conventions must be established to facilitate the com- 
parison and discussion which will follow. It will be agreed upon, in 
the language of this report, that the time origin in all the waveforms 
included here is at the farthest left hand side of the figures. Also, 
positive time is measured to the right hand side of the origin. The 
edge of a positive half -waveform closer to the time origin will be 
referred to as the "leading edge". On the other hand, the edge of a 
positive half-waveform farther from the time origin will be referred 
to as the "trailing edge". The waveform for various currents and 
voltages which are obtained from computer simulations will be referred 
to by the abbreviation CSWF. It should be pointed out that various 
CSWFs and their corresponding oscillogram counterparts have been ob- 
tained under slightly different operating conditions (speed, current 
command, etc.). Also it must be made clear that these investigators 
had no control over the experimental test procedures and conditions, 
since they were performed at a different facility. 
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5.1 VERIFICATION OF THE POWER CONDITIONER AND MACHINE MODEL 

The power conditioner and machine model is verified by compar- 
isons of simulation runs with actual oscillograms obtained from tests 
performed on the EMA. These comparisons are made for each of the 
three modes, as follows: 

Motoring 5000 rpm ICMD = 31.5 

Regeneration 5101 rpm ICMD = -38 

Plugging 240 rpm ICMD =32 

The input data used for the power conditioner and machine is given in 

Table (5.1-1). 

Computer simulations were obtained for the three basic modes of 
operation for the brushless DC machine-power conditioner package. The 
first mode is the forward motoring case at a rotor speed of 5000 r.p.m. 
The corresponding current and voltage waveforms through and across the 
numerous elements of the network for various commutation angles are 
given in Appendices A.l, A. 2, and A. 3, respectively. The second mode 
given here, is the generator (regenerative braking) case at a rotor 
speed of 5101 r.p.m. The corresponding current and voltage waveforms 
are given in Appendix B.l. The third mode studied is the forward 
braking case at a rotor speed of 240 r.p.m. The various waveforms 
obtained during this mode are displayed in Appendices C.l, C.2, and 
C.3. For the sake of convenience and comparison purposes some of the 
figures given in these appendices are repeated here with their corre- 
sponding experimentally obtained counterparts. 

Figure (5.1-1) shows experimentally obtained oscillograms of 
the machine phase currents at a speed of 5000 r.p.m. and with a comm- 




TABLE 5.1-1 The Values of the Parameters for the 
Power Conditioner and Machine Model 


R1 

.0333333/333.3333 (ON/OFF) 

R2 

.0333333/333.3333 

R3 

.0333333 

R4 

.0333333/333.3333 

R5 

.0333333/333.3333 

R6 

.0333333/333.3333 

R7 

.0596 

R8 

.0596 

R9 

.0596 

RIO 

.0333333 

Rll 

.0333333/333.3333 

R12 

.0333333/333.3333 

R13 

.0333333/333.3333 

R14 

.0333333/333.3333 

LI 

.000286 

L2 

.000486 (NONLINEAR) 

L3 

.000286 

L4 

.000286 

Cl 

.00099 

E4 

270 


UNITS: Ohms, Henries, Farads, and Volts. 



utation advance of 30 electrical degrees* The corresponding CSWF of 
the currents through phase A for commutation advance angles of 30, 22, 
and 38 electrical degrees are shown in Figures (5.1-3), (5.1-4), and 
(5,1-5). Figure (5.1-2) shows the spatial relationship between the 
output AA of the RPS for phase A and the command to inject current 
into phase A, Q (1) , for the chosen commutation angles. Examination 
of the current waveforms reveals the same characteristic over-shoot 
and under-shoot at the leading and trailing edges of the oscillogram 
and the CSWF. The positive peaks of the current oscillograms shown 
in Figure (5.1-1) are indicated by an (X) on the corresponding CSWF, 
Figure (5.1-3). Not:.cy the excellent agreement between the measured 
and computed current peaks. The transients in Figure (5.1-3) may 
appear exaggerated in comparison with Figure (5.1-1), but this is due 
to the "horizontal squeeze" and "vertical stretch" of the CSWF caused 
by the horizontal and vertical scaling of the plot. Both the oscill- 
ogram and the CSWF display the same chopper frequency of nine peaks 
per half cycle. Notice also the two slope current transient present 
in both cases in the middle of the positive and negative current cycle. 
Figures (5.1-4) and (5.1-5) demonstrate the effects of varying the 
commutation angle on the nature of the phase currents. This commutation 
angle also has a great effect on the electromagnetic torque of the 
machine as shown in Figures (A. 1-1), (A. 3-1), and (A. 2-1), respectively. 

The lower trace in Figure (5.1-6) depicts the oscillogram of 
the current through one of the transistor-diode switches in the six 
legged inverter-converter bridge leading to the machine phases. A 
corresponding CSWF is given in Figure (5.1-7). Excellent correlation 
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between the various ripples and pulses in the oscillogram and the CSWF 
are self evident upon examining these figures. Time congruence between 
the occurance of the various ripples and pulses is also evident in 
these figures. 

Figure (5.1-8) depicts the oscillogram of the main DC line 
current for a motoring mode. Notice the current dips which are occuring 
at regular intervals. The almost identical ripple and dip pattern is 
also evident in the CSWF given here in Figure (5.1-9). Again the 
agreement between the numerical and test results is excellent. 

Figure (5.1-10) gives the phase current for the generating (or 
regenerative braking) mode at a rotor speed of 5101 r.p.m., Figure 
(5.1-11) depicts the corresponding CSWF of one of the three phase 
currents at a constant rotor speed of 5101. Again agreement between 
the oscillogram and the CSWF is excellent. 

The last oscillogram, Figure (5.1-12) was made while the 
machine-power conditioner was operating in the braking mode at a 
speed of 240 r.p.m. The upper trace of this oscillogram shows a 
typical phase current for this mode of operation. The corresponding 
CSWF is given in Figure (5.1-13). The lower trace of Figure (5.1-12) 
shows the waveform of the current through one of the transistor switches 
connected to the permanent magnet machine. The corresponding CSWF 
is given in Figure (5.1-17). The CSWF shown in Figure (5.1-13) was 
obtained by using a reduced emf waveform, see Reference (15) . 

The CSWFs of Figures (5.1-14), (5.1-15), and (5.1-16) where obtained 
using the magnetized and demagnetized back emfs described earlier in 
Section 2.4. The reduced emfs are better in this case because 
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the large circulating phase currents result in the saturation of the 
stator teeth and yoke thereby increasing the ramf drop across the stator 
and hence reduced back emfs. The inclusion of this 

nonlinear stator reluctance in the back emf routine can be accomplished 
by a Finite Element analysis of the two dimensional nonlinear magnetic 
fields which govern the electromagnetic behavior of this machine. 

5.2 VERIFICATION OV THE TENTH AND FOURTEENTH ORDER EMA MODELS 

The overall EMA was modeled using both a detailed and a simpl- 
ified model of the power conditioner and machine network. Delco 
electronics. References (2), (3), and (4), performed a series of step 
responses on the actual hardware using four flap position commands 
of 1.1, 1.65, 2.2 and 2.75 flap degrees. The results of step responses 
are tabulated in table (5.2-1). This table contains the times it 
took the flap to reach the commanded flap positions for the first 
time. The maximum overshoot is also given. Close examination of this 
data reveals the excellent agreement between experimental and simu- 
lation results. 

The results of the 14th order EMA model are given in Appendix 
(D.l) . Figures (D.l-1) through (D.l-14) display some typical wave- 
forms produced by the EMA in response to a step position command. 

Figure (D.l-1) displays the status of the various EMA logic signals 
versus time. The extreme variations in rotor speed during such a 
maneuver is illustrated by the changing width of the RPS outputs 
(AA, BB, and CC) . Since each RPS pulse corresponds to 45 degrees of 
rotor displacement, the narrow RPS pulses occur at high rotor speeds 




and vice versa. The machine torque during the step response is given 
in Figure (D.l-2). The rapid fluctuations in the torque are due to 
the sawtooth component of the machine phase currents which is produced 
by chopper. Notice that as the machine velocity reaches its maximum 
(at approximately 0.15 sec) value, the average machine torque is signi- 
ficantly reduced. This is due to the increasing machine back emfs 
with increasing speed. This torque reduction at higher speeds results 
in velocity saturation which is graphically illustrated in Figure 
(5.2-1) (Figure (D.l-10)) and in Figure (5.2-2) which shows the corre- 
sponding observed velocity, Reference (4). 

The response of the flap to the step input is given in Figure 
(5.2-3). This waveform agrees with the oscillogram of the actual 
response given in Figure (5.2-4). Because of the linear relationship 
between machine current, IM„ and torque, TM, the machine current, 

Figure (D.l-4) is very similar to Figure (D.l-2). 

The computed velocity error given in Figure (D.l-5) shows 
how the velocity error approaches zero in time. The position error 
also appro* ches zero towards the end of the step response as expected, 
see Figure (D.l-14). 

The filtered current command IMC is displayed in Figure (D.l-6). 
The rectangular notch at t = .24 seconds is due to the lower current 
limit during plugging. The rate and magnitude limited current ICMD1/5 
is shown in Figure (5.2-5) as well as in Figure (D.l-7). This compares 
favorably with the actual current command waveform given in Figure 
(5.2-6). 

The action of the current rate limiter is displayed in Figure 


214 


(D.l-9). The reaction torque, TACT, given in Figure (D.l-11) shows the 
spring loading effect due to the nonzero value of the final flap 
position. 

In addition to these runs, the step response to the various 
position commands of Table (5.2-1) have been determined using a simpler 
10th order EMA model. Comparison of these results to test results 
demonstrates that the simpler model does an excellent job of predicting 
the step response. 

5.3 CONCLUSIONS AND RECOMMENDATIONS 

A comparison was made between the experimental test results. 
Reference (2) , and the numerical results presented in Appendices (A) 
through (E) . In all cases excellent agreement was found between the 
waveforms obtained from the actual machine-power conditioner package 
and those obtained from the mathematical (numerical) model. This 
close agreement between the behavior of the actual hardware and the 
mathematical model strongly suggests that the underlying assumptions 
and simplifications used in the development of this model are va2id 
and sound. 

Future improvements in this model can include the following 
endeavors: 

1. The development of a rigorous magnetic field analysis model 
for the determination of airgap flux density waveforms under load and 
no load. This model would constitute the heart of a program for the 
calculation of the machine back emfs in the machine-power conditioner 
network, as well as the calculation of the electromagnetic dynamic 
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torque. Such a package would allow the calculation of more accurate 
motor leakage inductances, which are necessary for accurate commutation 
studies. Such a model can also be used for the determination of better 
magnet-airgap configurations for optimum torque and dynamic characteristics. 

2. Increasing the order of the mathematical model in order to 
include the effects of the various capacitors and inductors which have 
been neglected in the present model. 

3. The incorporation of items (1) and (2) mentioned above into an 
overall EMA design optimization program. Such a program would determine 
the optimum system parameters given a set of system constraints. This 
type of capability would be useful in tailoring the system to a given 
application. 


APPENDIX A.l CALCOMP PLOTS: MOTORING 5000 rpm 30° ADVANCE 


RUN SUMMARY 

MODE: MOTORING 

PROGRAM MODE: PGMODE - 1 

CURRENT COMMAND: 1CMD =31.5 amps 

TIME STEP: TSNET = 5 ys 

•COMMUTATION SHIFT: SHIFT =0.0 mech rad 

CURRENT THRESHOLD: ITOL =4.0 amps 

ROTOR VELOCITY: (5000 rpm) (523.599 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 901 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW = 349 


TOTAL NUMBER OF INTEGRATIONS: NUMINT = 1800 
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Figure A. 1-16 
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APPENDIX A. 2 CALCOMP PLOTS: MOTORING 5000 rpm 38° ADVANCE 


RUN SUMMARY 

MODE: MOTORING 

PROGRAM MODE: PGMODE = 1 

CURRENT COMMAND: ICMD = 31.5 amps 

TIME STEP: TSNET = 5 ps 

COMMUTATION SHIFT: SHIFT = .034906 mech rad 

ITOL: 4.0 amps 

ROTOR VELOCITY: (5000 rpm) (523.599 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 901 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW = 359 


TOTAL NUMBER OF INTEGRATIONS: NUMINT = 1800 
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APPENDIX A. 3 CALCOMP PLOTS: MOTORING 5000 rpm 22° ADVANCE 


RUN SUMMARY 

MODE: MOTORING 
PROGRAM MODE: PGMOD = 1 

CURRENT COMMAND: ICMD =31.5 

TIME STEP: TSNET = 5 ps 

COMMUTATION SHIFT: SHIFT = -0.034906 mech rad 

CURRENT THRESHOLD: ITOL =4.0 amps 

ROTOR VELOCITY: (5000 rpm), (523,599 rad/sec) 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 1800 

NUMBER OF POINTS PLOTTED: NREC = 901 

TOTAL NUMBER OF NETWORK CHANGES: • NNETSW = 367 


EMfl LOGIC SIGNALS 
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Figure A. 3-1 
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APPENDIX B.l CALCOMP PLOTS: REGENERATION 5101 rpm 

RUN SUMMARY 

MODE: REGENERATION 

PROGRAM MODE: PGMODE = 1 

CURRENT COMMAND: ICMD = -38.0 amps 

TIME STEP: TSNET = 5 ys 

COMMUTATION SHIFT: SHIFT =0.0 mech rad 

CURRENT THRESHOLD: ITOL =4.0 amps 

ROTOR VELOCITY: (5101 rpm), (534.1756 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 901 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW = 415 


TOTAL NUMBER OF INTEGRATIONS: NUMINT = 1800 
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APPENDIX C.l CALCOMP PLOTS: PLUGGING 240 rpm 30° ADVANCE 


RUN SUMMARY 

MODE: PLUGGING 

PROGRAM MODE: PGMODE = 1 

CURRENT COMMAND: ICMD - 32.0 amps 

TIME STEP: TSNET = 10 ys 

COMMUTATION SHIFT: SHIFT =0.0 mech rad, 0.0 mech deg 

CURRENT THRESHOLD: ITOL =4.0 

ROTOR VELOCITY: (240 rpm), (24.13274 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 801 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW = 360 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 8000 


EMH LOGIC SIGNRLS 
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APPENDIX C.2 CALCOMP PLOTS: PLUGGING 240 rpm 38° ADVANCE 


RUN SUMMARY 

MODE: PLUGGING 

PROGRAM MODE: PGMODE = 1 

CURRENT COMMAND: ICMD = 32.0 amps 

TIME STEP: TSNET = 10 ps 

COMMUTATION SHIFT: SHIFT = 0.0349066 mech rad 

CURRENT THRESHOLD: ITOL =4.0 amps 

ROTOR VELOCITY: (240 rpm) , (24.3327 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 801 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW = 373 


a*,v 


TOTAL NUMBER OF INTEGRATIONS: NUMINT = 8000 






































APPENDIX C.3 CALCOMP PLOTS: PLUGGING 240 rpm 22° ADVANCE 


RUN SUMMARY 

MODE: PLUGGING 

PROGRAM MODE: PGMODE = 1 

CURRENT COMMAND: ICMD =32.0 amps 

TIME STEP: TSNET - 10 ys 

COMMUTATION SHIFT: SHIFT = -0.0349066 mech rad 

CURRENT THRESHOLD: ITOL = 4 

ROTOR VELOCITY: (240 rpm), (25.1327 rad/sec) 

TOTAL NUMBER OF POINTS PLOTTED: NREC = 801 

TOTAL NUMBER OF NETWORK CHANGES: NNETSW =372 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 8000 


EMfl LOGIC 5IGNRL5 



PDHI 
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APPENDIX D.l CALCOMP PLOTS: 14 th Order EMA Model DC = 2.75° 

RUN SUMMARY 

MODE: 14 th ORDER EMA MODEL 

PROGRAM MODE: PGMODE * 3 

NUMBER OF MACHINES; NMACH = 2 

DEFLECTION COMMAND: DC = 2.75 flap deg 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 40000 

TOTAL NUMBER OF MECHANICAL-POSITION LOOP INTEGRATIONS: NINTI = 2000 

TOTAL NUMBER OF VELOCITY LOOP INTEGRATIONS: 8000 

TOTAL SIMULATION TIME: 0.4 sec 

TOTAL NUMBER OF POINTS PLOTTED: 1001 

TOTAL NUMBER OF NETWORK CHANGES: 5811 

TIME TO DC: .179 sec 

TIME TO FIRST PEAK: .236 sec 

PEAK DEFLECTION: 3.233 flap deg 

MAXIMUM ROTOR SPEED: 10175 rpra 















12 





EMfl MOUNT REACTION TORQUE (TACT-XM(2S 



Figure D.l-11 



FLRP POSITION (XMm) 



1 f I ME IN SECONDS 
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APPENDIX E.l CALCOMP PLOTS: 10 th Order EMA Model DC - 1.1° 

RUN SUMMARY 

MODE: 10 th ORDER EMA MODEL 

PROGRAM MODE: PGMODE * 2 

NUMBER OF MACHINES: NMACH = 2 

DEFLECTION COMMAND: DC = 1.1 flap deg 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 2000 

TOTAL NUMBER OF MECHANICAL-POSITION LOOP INTEGRATIONS: NINTI - 2000 

TOTAL NUMBER OF VELOCITY LOOP INTEGRATIONS: 2000 

TOTAL SIMULATION TIME: 0.4 sec 

TOTAL NUMBER OF POINTS PLOTTED: 1001 

TOTAL NUMBER OF NETWORK CHANGES: 0 

TIME TO DC: .114 sec 

TIME TO FIRST PEAK: .163 sec 


PEAK DEFLECTION: 1.3329 flap deg 


EMH LOGIC SIGNRLS 



TIME IN SECONDS 



































CURRENT RATE LIMITER <XV(5» 
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APPENDIX E.2 CALCOMP PLOTS: 10 th Order EMA Model DC * 1.65° 

RUN SUMMARY 

MODE: 10 th ORDER EMA MODEL 

PROGRAM MODE: PGMODE = 2 

NUMBER OF MACHINES: NMACH = 2 

DEFLECTION COMMAND: DC = 1.65 flap deg 

TOTAL NUMBER OF INTEGRATIONS: NUMINT - 2000 

TOTAL NUMBER OF MECHANICAL-POSITION LOOP INTEGRATIONS: NINTI * 2000 

TOTAL NUMBER OF VELOCITY LOOP INTEGRATIONS: 2000 

TOTAL SIMULATION TIME: 0.4 sec 

TOTAL NUMBER OF POINTS PLOTTED: 1001 

TOTAL NUMBER OF NETWORK CHANGES: 0 

TIME TO DC: .139 sec 

TIME TO FIRST PEAK: .204 sec 

PEAK DEFLECTION: 2.14 flap deg 


EMfl LOGIC SIGNAL 5 



TIME IN SECONDS 



































Figure E.2-4 
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APPENDIX E.3 CALCOMP PLOTS: 10 th Order EMA Model DC - 2.2° 

PROGRAM SUMMARY 

MODE: 10 th ORDER EMA MODEL 

PROGRAM MODE: PGMODE = 2 

NUMBER OF MACHINES: NMACH = 2 

DELFECTION COMMAND: DC = 2.2 flap deg 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 2000 

TOTAL NUMBER OF MECHANICAL-POSITION LOOP INTEGRATIONS: NINTI - 2000 

TOTAL NUMBER OF VELOCITY LOOP INTEGRATIONS: 2000 

TOTAL SIMULATION TIME: 0.4 sec 

TOTAL NUMBER OF POINTS PLOTTED: 1001 

TOTAL NUMBER OF NETWORK CHANGES: 0 

TIME TO DC: .163 sec 

TIME TO FIRST PEAK: .229 sec 


PEAK DELFECTION: 2.73 flap deg 


EIMfl LOGIC SIGNALS 



IN SECONDS 
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TIME IN SECONOS 
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APPENDIX E. 4 CAT (COMP PLOTS: 10 th Order EMA Model DC ■ 2.75° 

RUN SUMMARY 

MODE: 10 th ORDER EMA MODEL 

PROGRAM MODE: PGMODE - 2 

NUMBER OF MACHINES: NMACH = 2 

DEFLECTION COMMAND: DC = 2.75 flap deg 

TOTAL NUMBER OF INTEGRATIONS: NUMINT = 2000 

TOTAL NUMBER OF MECHANICAL-POSITION LOOP INTEGRATIONS: NINITI = 

TOTAL NUMBER OF VELOCITY COOP INTEGRATIONS: 2000 

TOTAL SIMULATION TIME: 0.4 sec 

TOTAL NUMBER OF POINTS PLOTTED: 1001 

TOTAL NUMBER OF NETWORK CHANGES: 0 

TIME TO DC: .187 sec 

TIME TO FIRST PEAK: .253 sec 


PEAK DEFLECTION: 3.26 flap deg 


EMfl LOGIC SIGNALS 



TINE IN SECONDS 
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